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STABLE SIMPLEX SPLINE BASES 
FOR C 3 QUINTICS ON THE POWELL-SABIN 12-SPLIT 

TOM LYCHE AND GEORG MUNTINGH 


Abstract. For the space of C 3 quintics on the Powell-Sabin 12-split 
of a triangle, we determine explicitly the six symmetric simplex spline 
bases that reduce to a B-spline basis on each edge, have a positive parti¬ 
tion of unity, a Marsden identity that splits into real linear factors, and 
an intuitive domain mesh. The bases are stable in the L a0 norm with a 
condition number independent of the geometry, have a well-conditioned 
Lagrange interpolant at the domain points, and a quasi-interpolant with 
local approximation order 6. We show an h 2 bound for the distance be¬ 
tween the control points and the values of a spline at the corresponding 
domain points. For one of these bases we derive C°, C 1 , C 2 and C 3 con¬ 
ditions on the control points of two splines on adjacent macrotriangles. 


1. Introduction 


Piecewise polynomials or splines defined over triangulations form an indis¬ 
pensable tool in the sciences, with applications ranging from scattered data 
fitting to finding numerical solutions to partial differential equations. See 
[3,|l3j for comprehensive monographs. 

In applications like geometric modelling [4] and solving PDEs by isogeo¬ 


metric methods 10 one often desires a low degree spline with C 1 , C 2 or C 3 


smoothness. For a general triangulation, it was shown in Theorem l.(ii) of 
1261 that the minimal degree of a triangular C r element is 4r+l, e.g., degrees 
5, 9,13 for the classes C 1 , C 2 , C 3 . To obtain smooth splines of lower degree 
one can split each triangle in the triangulation into several subtriangles. One 
such split is the Powell-Sabin 12-split A of a triangle A; see Figure [ 2 ] On 

and 


18 


this split global C 1 smoothness can be obtained with degree only 2 
C 2 smoothness with degree only 5 112], 14,22 on any (planar) triangulation. 


Once a space is chosen one determines its dimension. The space 5^ and 
S 3 of C 1 quadratics and C 3 quintics on the 12-split of a single triangle 
have dimension 12 and 39, respectively. For a general triangulation T of a 
polygonal domain, we can replace each triangle in T by its 12-split to obtain 
a triangulation 7i2- The dimensions of the corresponding C 1 quadratic and 
C 2 quintic spaces (the latter with O 3 supersmoothness at the vertices of T 
and interior edges of the 12-split of each triangle in T) are 3|V| + \£\ and 
10|V| + 3|£|, where |V| and \£\ are the number of vertices and edges in T. 
Moreover, in addition to giving C 1 and C 2 spaces on any triangulation these 
spaces are suitable for multiresolution analysis |6l[7,14,17 . 


For a general smoothness, degree, and triangulation, it is a hard problem 
to find a basis of the corresponding spline space with all the usual properties 
of the univariate B-spline basis. One reason is that it is difficult to find a sin¬ 
gle recipe for the various valences and topologies of the triangulations. In 5 
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a basis, called the (quadratic) S-basis, was constructed for the C 1 quadratics 
on the Powell-Sabin 12-split on one triangle. The S-basis consists of sim¬ 
plex splines |15,191 and has all the usual properties of univariate B-splines, 
including a recurrence relation down to piecewise linear polynomials and 
a Marsden identity. Moreover, analogous to the Bernstein-Bezier case, C° 
and C 1 smoothness conditions were given tying the S-bases on neighboring 
triangles together to give C 1 smoothness on the refined triangulation 7i2- 
In the remainder of the next section we recall some background on splines 
and the 12-split. Section [3] introduces dimension formulas for spline spaces 
on the 12-split. In the next section some basic properties of simplex splines 
are recalled, after which an exhaustive list is derived of the C* 3 quintic sim¬ 
plex splines on the 12-split that reduce to a B-spline on the boundary. 
Section [5] introduces a barycentric form of the Marsden identity and de¬ 
scribes how the dual basis in [14 can be applied to find the simplex spline 


bases of 5| satisfying a Marsden identity. Making use of the computer alge¬ 
bra system Sage |20], these techniques are applied in Section [6] to discover 
six symmetric simplex spline bases that reduce to a B-spline basis on each 
boundary edge, have a positive partition of unity, a Marsden identity, and 
domain points with an intuitive control mesh and unisolvent for 5|. The 
concise and coordinate independent form of the barycentric Marsden iden¬ 
tity makes it possible to state these results in Table [3l Analogous to the 
Bernstein-Bezier case, we find C°, C 1 , C 2 and even C 3 conditions for one 
of these bases, tying its splines together on a general triangulation. One of 
the latter smoothness conditions involves just the control points in a sin¬ 
gle triangle and the geometry of the neighboring triangles, showing that C 3 
smoothness using 5 3 cannot be achieved on a general refined triangulation 
Ti 2 . Finally a conversion to the Hermite nodal basis from [14j is provided. 


2. Background 

2.1. Notation. On a triangulation T of a polygonal domain H C K 2 we 
define the spline spaces 

S r d {T) := {/ G <A(9) : /| A G n d VA G T}, r, d G Z, r > -1, d > 0, 

where 11^ is the space of bivariate polynomials of total degree at most d. 
With the convention 11^ := 0 and dimFO := 0 if d < 0, one has 

(1) dim 11^ = ^(d + 2)(d+ 1) + , d G Z, 

where z+ : = max{0, z} for any real number z. 

Any point x in a nondegenerate triangle [vi,v 2 ,v 3 \ can be represented 
by its barycentric coordinates (/?i, /3 2 > P 3 ), which are uniquely defined by 
x = P\V\ + ^ 2^2 + P 3 V 3 and /?i + P 2 + P 3 = 1. Similarly, each vector u is 
uniquely described by its directional coordinates , i.e., the triple — — 

f3' 2 , @3 — /3' 3 ) with (/3i,/ 3 2 ,/? 3 ) and (/3[, f3' 2 ,/3' 3 ) the barycentric coordinates of 
two points x and x' such that u = x — x'. Sometimes we write * as a linear 
combination of more than three vertices, in which case these coordinates are 
no longer unique. 
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Figure 1. Domain mesh for the quintic Bezier basis (left), 
the S-basis from J5] (middle) and the basis B c in this paper 
(right). 


A bivariate polynomial p of total degree d defined on a triangle Acl 2 
is conveniently represented by its Bezier form 

p ( x ) = c Bk B ijk( x ), B ijk ( x ) : = 

i+j-\-k=d 

i,j,k>0 


with (/?i, /? 2 , /3 3 ) the barycentric coordinates of x with respect to A. Here 
the Bfj k are the Bernstein basis polynomials of degree d and the coefficients 
Cijk are the Bezier ordinates of p. We associate each Bezier ordinate Cijk £ M 
to the domain point := ^1 + 3 ^ 2 + 2^3 £ M 2 and combine them into the 

control point c^k) £ M 3 . By connecting any two domain points 
and £i 2 j 2 k 2 by a line segment whenever \i\ — Z 2 1 + |ji — . 72 1 + |fci — fol = 1) one 
arrives at the domain mesh of p\ see Figure [laj The control mesh is defined 
similarly by connecting control points. 

We consider finite multisets K = {■u ”’ 1 • • • v™ n } C M 2 , in which the 
distinct elements v\,...,v n are counted with corresponding multiplicities 
mi,... ,m n > 0. Write |K| := m\ + ••• + m n for the total number of 
elements in K. For any two integers i,j, the Kronecker delta is the symbol 




1 if i = j, 

0 otherwise. 


For any set id 2 , define the indicator function 

U:R 2 — {0,1}. UW:={j|f Ha 


2.2. The Powell-Sabin 12-split. Given a triangle A = [v\,V 2 ,v^\ with 
vertices vi,v 2 ,v 3 write e\ := [^2,^3], e 2 := [v 3 , Vi], and e 3 := [^ 1 ,^ 2 ] 
for its (nonoriented) edges. Connecting vertices and the edge midpoints 
Vi := (v\ +i> 2 )/2, v$ := (t ?2 + ^ 3 )/2 and Vq := («i +u 3 )/2, we arrive at the 
Powell-Sabin 12-split A of A; see Figure [2a] for the labelling of the vertices 
V\,..., t?io and faces Ai,..., Ai 2 . 

To decide to which face of the 12-split points on the interior edges be¬ 
long, we follow the convention in 


which can be 


shown in Figure 2b 
quickly computed by Algorithm 1.1 in [5|. If K is a multiset satisfying 
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Figure 2. The Powell-Sabin 12-split with labelling of ver¬ 
tices and faces (left), and a scheme assigning every point in 
the macrotriangle to a unique face of the 12-split (right). 


K C {iq,..., Dio} as sets, its convex hull [K] is a union of some of the faces 
of the 12-split, and we define the half-open convex hull [K) as the union of 
the corresponding half-open faces depicted in Figure [2b| 

Remark 1. While in 126 it was shown that on a general triangulation de¬ 
gree 9 is needed to achieve global C 2 smoothness, let us show that on the 12- 
split T 12 of a triangulation the degree is necessarily at least 5. For r = 0,1, 2, 
a C r spline of degree r+2 on the knot multiset {0 r+3 0.5 2 l r + 3 } is determined 
by 5 + r degrees of freedom. To aim for C 2 smoothness using quartics, one 
uses the minimal number of degrees of freedom as follows. At each vertex 
fix derivatives up to order 2, and on each edge fix one additional value, two 
cross-boundary derivatives, and three second-order cross-boundary deriva¬ 
tives. Thus in total 36 degrees of freedom are needed, while the space 5 2 (A) 
only has dimension 34; see Table [l] 

2.3. A basis for the dual space of a space of C 3 quintics. Let A 

be the 12-split of a triangle A with vertices V = {di,D 2 ,d 3 } and edges 
8 = {[di,d 2 ], [d 2 ,d 3 ], [d 3 ,di]}. For any edge e = [vi,Vj\ G 8 with opposing 
vertex v let 

3 Vi + Vj Vi + Vj Vi + 3d 7 - 

Ql,e ■= - m e ■= - 2 “^ 92,e : = - 

be its midpoint and quarterpoints. With e v the point evaluation at v and 
D u the directional derivative in the direction u, let 


( 2 ) 


A *— U 




D 2 F 

Ql 


Du e ,£q 2 D u }'i 


i+j <3 eG£ 

i,j> 0 


see Figure [3} Here, for every vertex v and edge e, the symbols x v ,y v are 
any linearly independent vectors and u e is any vector not tangent to e, 
for example the outside unit normal as shown in Figure [3} It was shown 
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Figure 3. A basis for the dual space of on a single tri¬ 
angle. A bullet represents a point evaluation, three circles 
represent all derivatives up to order three, and a single and 
double arrow represent a first- and second-order directional 
derivative. These derivatives are evaluated at the rear end of 
the arrows, which are located at the midpoints and quarter- 
points. 


m 


14, Theorem 4] that A is a basis for the dual space to 5|(A), which 


therefore has dimension 39. 


3. Dimension formulas 


Consider a polygonal domain !lci 2 with a triangulation T with sets of 
vertices V, edges £, faces T, and 12-split refinement 7i2- Let 


S 5 2 ’ 3 (Ti 2 ) := {s € 5|(Ti 2 ) : s G C 3 (A)nC 3 (u) for all A € A, v G V}. 

Here s G C 3 (v ) means that all polynomials such that A is a triangle 
with vertex at v have common derivatives up to order three at the point v. 
Note that if T consists of a single triangle, then 5 2 ’ 3 (7i2) = 5 3 (7i2)- 

Since A from © specifies the value and partial derivatives up to or¬ 
der three at each vertex of T and the value, first- and second-order cross¬ 
boundary derivatives at each edge of T, the following theorem is an imme¬ 
diate consequence of 14, Theorem 4], 


Theorem 2. For any triangulation T with |V| vertices and \£\ edges, the 
set A is a basis for the dual space of S^’ (7i2)- Ln particular 


dimS 2 ’ 3 (Ti 2 ) = 10|V| + 3|£|. 


Next, let A denote the 12-split triangulation of a single triangle A. For 
future reference, we state the following formula for the dimension of the 
space of C r splines of degree d on A, which is a special case of Theorem 3.1 
in |2| (also compare 12l]). 
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dim Sj (A) 

c - 1 

C° 

C 1 

C 2 

C 3 

C 4 

C 5 

C 6 

C 7 

C 8 

C 9 

d = 0 

12 

1 










d= 1 

36 

10 

3 









d = 2 

72 

31 

12 

6 








d = 3 

120 

64 

30 

16 

10 







d = 4 

180 

109 

60 

34 

21 

15 






d = 5 

252 

166 

102 

61 

39 

27 

21 





d = 6 

336 

235 

156 

100 

66 

46 

34 

28 




d = 7 

432 

316 

222 

151 

102 

73 

54 

42 

36 



d = 8 

540 

409 

300 

214 

150 

109 

81 

63 

51 

45 


d = 9 

660 

514 

390 

289 

210 

154 

117 

91 

73 

61 

55 


Table 1. Dimensions of 5^(A), with (r, d) = (3,5) highlighted. 


Theorem 3. For any integers d,r with d> 0 and d> r > —1, 
dim S r d{ 


= ^(r + l)(r + 2) + ^(d — r)(d — r + 1) 


(3) 


d—r 


+ -(d-2r~ l)(d - 2r)+ + ^(r - 2j + 1)+. 


j'=i 


To quickly look up dim 5^ (A) for small values of r and d, we have listed 
these first dimensions in Tabled] 


4. Simplex splines 


In this section we first recall the definition and some basic properties of 
the simplex spline, and then proceed to determine the C 8 quintic simplex 
splines on the 12-split that reduce to a B-spline on the boundary. For a 
comprehensive account of the theory of simplex splines, see 


15,19 


4.1. Definition and properties. The following definition of the simplex 
spline is convenient for our purposes. 


Definition 1. For any finite multiset K = {v™ 1 ■ ■ ■ v™ 0 10 } C M 2 composed 
of vertices of A, the (area normalized) simplex spline Q[ K] : M 2 —> M is 
recursively defined by 


Q[ K](*) 


0 if area([K]) = 0, 

< 1 [K)( a; ) ^ a ([K]) if area([K]) ^ 0 and |K| = 3, 

T,]liPjQ[ K \ v j\( x ) if area([K]) / 0 and |K| > 3, 


with x = /3\Vi + ■ ■ - + /3ior>io, Pi + ■ ■ ■ + /3io = b and (3i = 0 whenever rrq = 0. 


By Theorem 4 in Jl5| this definition is independent of the choice of the 
j3j. It is well known that Q[K] is a piecewise polynomial with support 
[K] C A and of total degree at most |K| — 3. One shows by induction 

on |K| that f R2 Q\K]{x)dx = area(A) • (^ K i _1 ) . Although M[ K] := 

area(A) -1 (^ k ^~ 1 )<5[K] has unit integral and is used more frequently, our 
discussion is simpler in terms of Q [K]. 
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Whenever m? = mg = mg = m±g = 0, we use the graphical notation 

:= QlvivivWtvfv®. 

Example 1. The linear S-spline basis in [ 5 ] only uses vertices tq ,... ,vq, 
V\q. while the quadratic S-spline basis only uses V\,... ,Vq. It is given by 


(4) 


Sj, 2 = 


area([Kj j2 ]) 


M[ K ii2 ] = 


area([Ky 2 ]) 


6 L J '* i area(A) 

where by (2.5) in |5|, as (unordered) sets, 

{g[K 1 ) 2 ],...,Q[K 12)2 ]} = { 

(5) 


Q[ Kd, j = 1, - - ■, 12, 


}■ 


Example 2. If K = Vj j+1 v^ k+1 } has three distinct elements with 

area([K]) > 0, then, with (Pi,Pj,Pk) the barycentric coordinates of x with 
respect to [vi,Vj,Vk], it follows by induction that 


Q[ K](a) 


4[K)(®) 


area ( A ) + H + ^Y- m 8ftk 

area([K]) m\nj\fj, k \ 1 j k 


is, up to a scalar, a Bernstein polynomial on [K). 


Continuity. For any edge e of A, if K has at most m knots (counting multi¬ 
plicities) along the affine hull of e, then Q[K] is |K|— m —2 times continuously 
differentiable over e. For instance, every C 3 quintic simplex spline on A has 
at most three knots along the affine hull of any interior edge e. 


Differentiation. Let K = { vf 1 ■ ■ ■ vf^ 0 } be a finite multiset. If u = aq-iq + 
■ ■ ■ + aio^io is such that aq + • • • + aqo = 0 and oq = 0 whenever mj = 0, 
then one has a differentiation formula 

10 

(6) D U Q[ K] = (|K| - 3) £ aj Q[K\ Vj }. 

3 =1 

Knot insertion. If y = f3\V\ + • • • + f3igV\g is such that (3 1 + • • • + f3 10 = 1 
and flj = 0 whenever mj = 0, then one has a knot insertion formula 

10 

(7) Q[K} = Y,PjQ{(KUy)\vj]. 

3 =1 

For instance, if Vi, v 2 6 K, then, since V 4 = \v\ + |v 2 , 

(8) Q[ K] = l -Q[(K U « 4 )\vi] + \Q[{ K U v 4 )\t> 2 ], 

and for example 
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Restriction to an edge. Let e = [vi,Vk] be an edge of A with midpoint Vj 
and let <pik(t) ■= (1 — t)vi + tv k ■ By induction on |K|, 


(9) Q[K]otp ik (t) 


0 


area(A) 
area( [K]) 


B(t) 


if rrii + nij + nik < |K| — 1, 
if m,; + rrij + m k = |K| — 1, 


where B is the univariate B-spline with knot multiset {0 mi 0.5 mj l mfe }. 

We say that Q[ K] reduces to a B-spline on the boundary when B is one 
of the consecutive univariate quintic B-splines B\,,B% on the open knot 
multiset {0 6 0.5 2 l 6 }; see Table [dj Similarly a basis B = {Si,..., S 39 } of 
Sf (A) reduces to a B-spline basis on the boundary when 

{Sio^fc, ..., S 39 oip ik ] = {(Bf) 1 ■ ■ ■ (Bf ) 1 0 31 }, 1 < i < k < 3, 

as multisets. This scaling of B ensures simple C° conditions for connecting 
two adjacent patches expressed in terms of B. 


Symmetries. The dihedral group S 3 of the equilateral triangle consists of 
the identity, two rotations and three reflections, i.e., 





The affine bijection sending v k to (cos 27r/c/3, sin 2nk/3), for k = 1,2,3, 
maps A to the 12-split of an equilateral triangle. Through this correspon¬ 
dence, the dihedral group permutes the vertices Vi,...,Vio of A. Every 
such permutation a induces a bijection Qlv™ 1 ■ ■ ■ v^ 10 ] 1 —> Q[cr(ui) mi • • • 
o-(fio) mi °] 011 the set of all simplex splines on A. For any set B of simplex 
splines, we write 

[B]s 3 := |Qk(K)] : Q[ K] G B, a G S 3 } 

for the S 3 equivalence class of B, i.e., the set of simplex splines related to B 
by a symmetry in S 3 . For example, 



and ([5]) takes the compact form 



We say that B is S^-invariant whenever [B]s 3 = B. 

4.2. C 3 quintic simplex splines on the 12-split. Any simplex spline 
Q[K] of degree d = 5 on A is specified by a multiset K = { vf 11 ■ ■ ■ "C^ 10 } 
satisfying 

( 10 ) mi + 7712 + • • • + 77710 = d + 3 = 8 . 
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Lemma 1. Suppose a quintic simplex spline Q[K] on A is of class C 3 . 
Then 

(11) mt = mg = m 9 = m io = 0, 
and 

(12) mi + 7715, 7713 + 7714, 7712 + 7776, 7774 + 7716 , 7714 + TO5, 7715 + 7716 < 3, 
whenever both multiplicities are nonzero. 


Proof. In the 12-split, certain knot lines do not appear, leading to the con¬ 
ditions 

771i?7lg = 771i?7lg = mg?71g = 0, 

(13) 77127717 = 7712 771g = 77l7?71g = 0, 

77 I 377 I 7 = ?7l3771g = 77l7?71g = 0. 

To achieve C* 3 smoothness over the remaining knot lines, the sum of the 
multiplicities along each line in the 12-split has to be at most three, 

7711 + 7715 + 7717 + 77110 < 3, 7714 + 7716 + 7717 < 3, 

(14) 7712 + 777-6 + 77lg + 77110 < 3, 777-4 + ?7l 5 + 771g < 3, 

7713 + 777-4 + 771g + ?7lio < 3, 777-5 + 7716 + 77lg < 3, 


whenever the line contains at least two knots with positive multiplicities. 
Suppose 7717 > 1. Then (13) implies 777-2 = 7713 = mg = mg = 0, and by 


(10) and (14), 8 = (mi + ms + m .7 + mio) + (m .4 + me) <3 + 2, which is a 
contradiction. It follows that 7717 (and similarly ms and mg) must be equal 
to zero. Moreover, by (10) and (14), 8 = (mi + 7715 + 7717 + m-io) + (m 2 + 
?7i6 + ms + mio) + (m 3 + 7714 + mg + 7?iio) — 2mio < 9 — 2mio. Therefore 
? 7 iio = 0 and (11) holds, and (|12|) follows immediately from (14). □ 


In addition we demand that Q[K] reduces to a B-spline on the boundary. 
By (§, if mi + mj + 77ifc < 7, then Q[K]| e = 0 and this condition is satisfied. 


The 

remaining case m 

■i + mj 

+ m k = 

: 7 yields the 

conditions 





not (mi 

+ 7714 

+ m 2 = 

= 7 and 

777,4 

> 

3), 








not (mi 

+ 7714 

+ m 2 = 

= 7 and 

mi 

> 

1 and 

m 2 

> 

1 and 

m .4 


2) 

/I C \ 

not (m 2 

+ rn 5 

+ m 3 = 

= 7 and 

m 5 

> 

3), 







(15) 

not (m 2 

+ rn 5 

+ m 3 = 

= 7 and 

m 2 

> 

1 and 

7713 

> 

1 and 

m- 5 


2) 


not (mi 

+ me 

+ m 3 = 

= 7 and 

m 6 

> 

3), 








not (mi 

+ m 6 

+ m 3 = 

= 7 and 

mi 

> 

1 and 

m 3 

> 

1 and 

m-6 


2) 


Adding the entries in the second column of (14) and using © gives 
2?7i4 + 2?7i5 + 2 me < 9, implying, by (10), 


( 16 ) mi + m2 + m3 > 4 , 7714 + 7715 + 77 i 6 < 4 . 


Theorem 4. ITif/i one representative for each S 3 equivalence class, Tabled 
is an exhaustive list of the C 3 quintic simplex splines on A that reduce to a 
B-spline on the boundary. 
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Q a 

Q b 

Q c 

Q d 

Q e 

Q f 

Q 9 

Q h 

Q i 

Q j 

& 

& 

A 

0%© 

©°©© 

^© 

©©© 



©1© 

Q k 

Q l 

Q m 

Q n 

Q° 

Q p 

Q q 

Q r 

Q s 

Q t 

©©© 


o©°© 

©s© 

©§© 

0©© 

©©© 

©©°© 

©^© 

A 


Table 2 . The C 3 quintic simplex splines on A, one rep¬ 
resentative for each S3 equivalence class, that reduce to a 
B-spline on the boundary. 


Proof. Recall from that 7717 = ms = 1^9 = Tfi 10 = 0 - We distinguish 
cases according to the support [K] of Q[K], up to a symmetry in S3. 

Case 0 , no corner included, [K] = [774,715, Tie]: By ( 10 ), 7774+777,5+7776 = 8, 
so that the sum of two of these multiplicities will be at least 5 , contradict¬ 


ing ( 12 ). Therefore this case does not occur. 


Case la, 1 corner included, [K] = [774,7)4,776]: For a positive support 
7774,7774, 7776 > 1 , and since 777,4 + < 3 by ([ 12 ]), we obtain 


[Q a ]s 3 = 


[</]s 3 = 


-I S3 


-I S3 


Case lb, 1 corner included, [K] = [774,774,775,7)6]: By ( 10 ) and ( 16 ) one 
has 7774 = 8 — 7774 — 7775 — 7776 > 4 , contradicting 7774 + 7775 < 3 from ©• 
Therefore this case does not occur. 

Case 2 a, 2 corners included, [K] = [vi,V2,ve]'. If 7774 G { 0 , 1 }, then 


777 6 > 2 by the second line in ( 15 ), and since mg < 3 — m2 by ( 12 ), it follows 


that m2 = 1 and 777 6 = 2 . We obtain 


= 


[Q d ]s 3 = 


J s 3 


s 3 


If 7774 > 2 , then, since 7774 + me < 3 by ( 12 ), one has 7774 = 2 and me = 1 . 
Since m2 < 3 — me = 2 by © , we obtain 


[Q e ]s 3 = 


-I S 3 


7 [Q J ]s 3 = 


Js 3 


Case 2 b, 2 corners included, [K] = [7)4,7)2,775,7)6]: Since 7775 + 7776 < 3 by 
( 12 ), one has 1 < 7775,7776 < 2 . Suppose me = 2 . Then m2 = 7775 = 1 and 


7774 < 1 by ( 12 ), and 7774 + 7775 = 8 — m2 — 7774 — me > 4 , contradicting ( 12 ). 


We conclude that me = 1 and similarly that 7775 = 1 . Then 7774 + 7772 + 7774 = 


8 — me — me = 6 , and since ( 12 ) implies 7774,777,2, m 4 < 2 , one obtains 


[Q 9 ]s 3 = 


J S 3 


Case 3 , 3 corners included, [K] = [7)4,772,773]: We distinguish cases for 
(7774,7775,7776), with 7774 > me > me, first by 7774, and then by 7774 + 7775 + me, 


which is at most 4 by ( 16 ). 
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( 0 , 0 , 0 ) 

One has m \, m 2 

m 3 


[Q h }s 3 = 


? 

S3 

( 1 , 0 , 0 ) 

One has 5 

< m 3 

< 3 


[Q J ]sa = 

cJL 

5 

S3 

( 1 , 1 , 0 ) 

One has mi, m 3 

< 2 


[Q Z ]s 3 = 

©€&> 

5 

S3 

( 1 , 1 , 1 ) 

One has mi, m 2 

m 3 


[Q 0 ]s 3 = 


S3 

( 2 , 0 , 0 ) 

One has m 3 = 1 

by ( 


[Q P ]Sa = 

©5© 

s 3 

( 2 , 1 , 0 ) 

One has m 3 = 1 

and 


[Q P ]s s = 

C*© 1 ) 

S3 

( 2 , 1 , 1 ) 

One has m 3 = 1 

and 


[Q*]s 3 = 


S 3 


, [Q l ]s 3 = 

777,4 by i 

, [Q k ]s 3 = 


- S3 


-I S3 


, [<T]s 3 = 


, [Q n ]s 3 = 


Sa 


> Kfls. = 


, [<3 S ]s 3 = 


JS 3 


JS 3 


-I S 3 


□ 


5. Simplex spline bases for Sf 

Let A as in ©> be a basis of the dual space of <S|(A). In this section we 
describe a recipe for determining the S^-invariant simplex spline bases that 
reduce to a B-spline basis on the boundary, having a positive partition of 
unity and a Marsden identity with only real linear factors. 

5.1. Potential bases. Of the C 3 quintic simplex splines in Table [ 2 J only 
Q a Q b Q c Q e Q f Q p Q q 



are nonzero on [ui ,v?\. Any S' 3 -invariant basis reducing to a B-spline basis 
on the boundary should therefore contain the S 3 equivalence class of one of 
the eight possible combinations in the Cartesian product 



This determines 3 + 6 + 6 + 6 = 21 of the 39 basis elements. 

Of the 13 remaining S 3 equivalence classes in Table [ 2 ] of simplex splines 
that are zero on [vi, V 2 ], there are 6 of size 3, 
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Q 9 Q h Q i Q l Q n Q° 



and 7 of size 6 , 


Q d Q j Q k Q m Q r Q s Q t 



For each of the above 8 choices of 21 simplex splines that are nontrivial on 
the boundary <9A, we complete the basis by adding S 3 equivalence classes 
with together 18 elements, resulting in 



potential Ss-invariant simplex spline bases B of 5 | that reduce to a El- 
spline basis on the boundary. One selects the linearly independent sets B 
by checking that the collocation-like matrix {A(Q)} ( g e g j AeA has full rank. 

5 . 2 . Positive partition of unity. Let B = {Qi, . . ., Q39} be an ordered 
simplex spline basis of (A) that reduces to a B-spline basis on the bound¬ 
ary. We desire that B has a positive partition of unity, i.e., that there exist 
weights w \,..., W39 > 0 for which 
39 

(17) y ^WjQj(x) = 1 

i=l 

holds identically. Applying the functionals in A = {Ai,..., A39} yields 
39 

^WiXjiQi) = Xj{l), j = 1,... ,39, 

i =1 

which has a unique solution (w 1,..., W39) £ M 39 by the linear independence 
of B and A. One then checks whether the weights Wi are all positive. 


5.3. Marsden identity. More generally, we would like B to satisfy a Mars- 
den identity 

39 

(18) (1 + x T y) 5 = ^2wiQi(x)tpi(y), x £ A, y £ M 2 , 

i= 1 

for certain dual polynomials 'ipi and dual points p* r of the form 


My) ■= n^+PiJs/), p tr e 


r = 1,..., 5, i = 1,..., 39. 


r =1 


In particular one recovers the partition of unity by setting y = 0. Similarly 
one can generate all other polynomials of degree at most 5. For instance, 
differentiating (18) with respect to y and evaluating at y = 0 yields 
39 


(19) x = ^T€ i w i Qi(x), 


1=1 


& := 


ph + ---+pu 


y =0 


where $ >i is the domain point associated to Qi. 
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While the Marsden identity (18) is the form commonly encountered in the 
literature [8j[9], we instead present a barycentric form that is independent 
of the vertices of the macrotriangle. 


Theorem 5 (Barycentric Marsden identity). Let f3j = f3j(x),j = 1,2,3, be 
the barycentric coordinates of x 6l 2 with respect to A= [v\,V 2 ,v^]. Then 


(18) is equivalent to 


39 

(20) (fiici + /3 2 c 2 + /3 3 c 3 ) 5 = 'y]w i Q i (PiVi + /3 2 v 2 + /0 3 v 3 )^i(ci, c 2 , c 3 ), 

i =1 

where x £ A, ci, C 2 , c 3 £ M, and, for i = 1,..., 39, 

5 

^i(ci,c 2 ,c 3 ) := (/3i(p* r )ci + /3 2 (p- ir )c 2 + /3 3 (p* r )c 3 ). 

r=l 


Proof. Let Xj := (1 ,Vj) T £ M 3 , j = 1,2,3. Since {Xi,X 2 ,X 3 } is linearly 
independent, there are Y ? ; £ M 3 , i = 1,2,3, such that XjY* = Sij. Given 
x £ A and ci,c 2 ,c 3 £ M we define X,Y £ M 3 by X := (1,*) T and Y = 
(1, y) T := ciYi + c 2 Y 2 + c 3 Y 3 . Since X = /3 1 (x)X 1 + /3 2 (*)X 2 + P 3 (x)X 3 , 


1 + a? T i/ = X t Y = (5i(x)c\ + /3 2 (*)c 2 + p 3 {x)c 3 , 
and in particular with a; = p* r and P* r := (l,p* r ) T , 

1 + Pijy = P*,r Y = Pi (p* r )c\ + (3 2 (p* r )c 2 + /3 3 (Pi,r) C 3- 


It follows that (20) is equivalent to (18). 


□ 


For a compact representation of the dual polynomials, we introduce the 
shorthands (compare Figure [2a|) 


( 21 ) 


c 4 := 
c 7 := 


ci + c 2 
2 : 
C4 + C6 


C5 := 

c 8 := 


c 2 + c 3 

2 ' 

C 4 + C 5 


c 6 : = 
c 9 : = 


ci + c 3 
2 ! 
C 5 + c 6 


cio : = 


ci + c 2 + c 3 


Example 3. The barycentric Marsden identity for the quadratic S-basis in 
5, Theorem 3.1] (cf. Example [I]) is 


(ci/3i + C 2 P 2 + c 3 /3 3 ) 2 = ci 


1 2 
+ 4 Cl 

+ 4 ci 

1 2 

+ “C 3 

4 d 


+ -C4C10 


+ ^C 5 Cio 


+ |C 6 Cio 


+ C 2 


+ 2 ClC4 


:c 2 c 5 


+ 2 C 3 C 6 


+ c 3 


+ 2 C2C4 , 


+ 2 C3C5 


+ 2 C 1 C 6 


where the weigths follow from Q. 
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For every basis B = {Qi,..., Q 39 } with a positive partition of unity 
we apply the functionals in A to (18) with respect to x, which gives 


/Ai(Qi) 

\A39(Qi) 


Ai (Q39) 


^39(^39) / 




( Ai ((1 + x T y) 5 ) 


\w 39 ip 3 9 (y) J \A 39 ((1 + x T y) 5 ) / 


This system has a unique solution in the module M[y] 39 , with each compo¬ 
nent vjiil>i(y) a polynomial of degree at most 5 by Cramer’s rule. The basis 


B has a Marsden identity (18) if and only if these polynomials split into real 
linear factors. 


6. Main results 


Some of the remaining computations are too large to carry out by hand. We 
have therefore implemented the above computations in Sage 20 . From the 


website [16] of the second author, the resulting worksheet can be downloaded 
and tried out online in SageMathCloud. 


6.1. Six bases for S |. Checking linear independence for each of the 3648 
potential bases, we discover that there are 1024 S 3 -invariant simplex spline 
bases that reduce to a B-spline basis on the boundary. There are 243 such 
bases with a nonnegative partition of unity, of which there are 47 bases with 
a positive partition of unity. Of these there are 9 with all domain points 
inside the macrotriangle, of which there are 7 with precisely 8 domain points 
on each boundary edge. Of these there are 6 bases B ai Bb,B c ,Bd,B e ,Bf 
for which each dual polynomial only has real linear factors. These bases, 
together with their weights dual polynomials T*, and domain points 
are listed in Table [3j For instance, the highlighted rows in the table yield 
the set 



We summarize these results in the following theorem. 


Theorem 6. The sets B = B a ,Bb, B c , Bd, B e , Bf are the only sets satisfying: 

(1) B is a basis o/5f(A) consisting of simplex splines. 

(2) B is S 3 -invariant. 

(3) B reduces to a B-spline basis on the boundary. 

(4) B has a positive partition of unity and a Marsden identity pp| ; for 
which the dual polynomials have only real linear factors; see Tabled 

(5) B has all its domain points inside the macrotriangle A, with precisely 
8 domain points on each edge of A. 


Remark 7. Table [3] shows that, for the simplex splines 


[Q 9 ,Q S ,Q% 3 = 



the weights and dual polynomials depend on the entire basis, and they 
cannot be determined directly from the corresponding knot multiset. 
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Q a 

Q b 

Q e 

Q f 

Q 9 

Q l 




& 

o°©o 

^© 

©©© 

©o°o 

B a 


l r 5 

2 cfc 4 

5 C 1 C 4 

5 C 1 C 2 C 4 

0 

C 2 c 4 c 5 



(1,0,0) 


(f.i.o) 

(f’!’ 0 ) 

X 

(- - -) 

U0> 5’ 10/ 

B b 


l r 5 

1 f A c . 
4 c 1 c 4 

| c l c 4 

(f.J.o) 

2 C 1 C 2 C 4 

0 

C 3 C 4 C5 



(1,0,0) 


(i>§>°) 

X 

f -L 1 JO 

U0> 5’ 10/ 

Be 


1 ,,5 

4 C 1 

l r 4 r . 

4 C 1 C 4 

1 „3,,2 
2 L 1 C 4 

(* 1 o') 

12 2 
2 C I C 2 C I 

|cic 2 c|cio 

c|c 4 c 5 



(1,0,0) 

( 9 1 o) 

( 3 2 o') 

( 7 7 1 ) 

( 1 1 M 


Vio > 10 > u / 

■ HHHHI - 

■ 4 5’ 1 H 

U5’ 15’ 15/ 

VlO’ 51 10/ 

B d 

uh’F* 

I r 5 

4 C 1 

4 C 1 C 4 

| c l c 4 

(U°) 

0 

0 

C 2 c 4 c 5 



(1,0,0) 


X 

X 

(-L 1 

U0> 5’ 10/ 

Be 

uh’F* 

I r 5 

4 C 1 

l r 4 rA 

4 C 1 C 4 

f c l c 4 

(M.°) 

0 

0 

C 2 c 4 c 5 



(1,0,0) 

(&4.o) 

X 

X 

(J_ 4 J_\ 
U0> 5’ 10/ 

Bf 


l r 5 

4 C 1 

1~4„ 

4 C 1 C 4 

§ c l c 4 

0 

iCic 2 c 3 c| 

C 2 c 4 c 5 



(1,0,0) 


(fi.o) 

X 

(2 2 14 

45’ 5’ 5/ 

(ih’S’ll)) 



Q n 

Q° 

Q q 

Q r 

Q s 

Q l 


© 0 © ©©© ©^^) ©^© 


Ba 

lO'F* 

|cic|c 3 cio 

0 


0 

^cic|c 4 c 5 

|cic|c 4 cio 

0 



1 4 

7 4 4 

X 


X 

( 3 3 1 4 

/n 17 14 

X 


4 15 ’ 

15’ 15/ 


110’ 5’ 10/ 

4 30’ 30’ 15/ 

B b 

114 'F* 


0 

|cic 2 c 3 c 4 cio 

0 

icic|c 4 c 5 

|cic|c 4 ci 0 

0 




X 

/II 11 

4 i 

X 

( 3 3 14 

/II 17 14 

X 



1.30 ’ 30’ 

15/ 

U0’ 5’ 10/ 

4 30 ’ 30 ’ 15/ 

Be 

Hh'F* 


0 

0 


0 

^cic|c 4 c 5 

0 

|ciC 2 C 4 C5Cio 




X 

X 


X 

(A. 3 JO 
U0 ’ 5 ’ 10 J 

X 

(U JL l) 

V 30 ’ 15 ’ 6 / 

B d 

114'F* 

|cic|c 3 cio 

0 


c 2 c 2 c| 

icic|c 4 C 5 

4ciC 2 C 3 C 4 

0 



/ 4 

7 4 4 

X 


(3 2 n 4 

( 3 3 1 4 

(3 114 

X 


V 15 5 

15’ 15/ 


(5’5’ U ) 

U0’ 5’ 10/ 

U0’ 2’ 5/ 

Be 

114 'F* 


0 

|cic 2 c 3 c 4 cio 

c 2 c 2 c| 

icic|c 4 C 5 

2 C 1 C 2 C 3 C 4 

0 




X 

(11 11 

4 i 

(3 2 n 4 

/ 3 3 1 4 

(3 114 

X 



4 30 ’ 30’ 

15/ 

U’ 5’ U ) 

V 10’ 5’ 10/ 

410’ 2’ 5/ 

Bf 

114 'F* 


0 

0 


c\c 2 c\ 

|cic|c 4 c 5 

0 

2 c ic 2 c 3 c 4 cs 




X 

X 


( 3 2 0) 

t 3 3 1 ) 

X 

( 7 2 14 




\5 ’ 5 ’ U ) 

V 10 ’ 5’ 10/ 

420’ 5’ 4/ 


Table 3 . For the six bases B a ,...,Bf and for each contained 
representative Q* in Table [2j the table lists the weight re* 
times the dual polynomial ’ll*, using the shorthands in (21), 
and the barycentric coordinates of the domain point 


Dual points and domain points. One immediately reads off the dual points 
from the dual polynomials in Table [3j simply by replacing ‘c’ by l v\ For 
instance, in each basis the simplex spline 
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Figure 4. For the six bases in Table [ 3 J the figure shows the 
domain points © and, for each S 3 equivalence class, the label 
of a representative in Table [2j 


has dual polynomial = C 2 C 4 C 5 , and therefore dual points {v^v^v^}. By 
(19) the corresponding domain point is ( 3 v 2 + V 4 + v 3 )/5 = + |i >2 + 

Thus one obtains for each basis all domain points, which are listed 
in Table [3] and shown in Figure |4j The domain points of the basis B c are 
connected to form the domain mesh in Figure [Tc} To preserve the symmetry 
of A, the domain points are forced to form a hybrid mesh with triangles, 
quadrilaterals, and a hexagon in the center. 


Polynomial reproduction. Following [11], we define for any nonnegative in¬ 
tegers ii,* 2,*3 the “coefficient of” operator 


C 1 c 2 c 3 


f=T 74 


d 5 f 


dcYdc^dcl 3 


( 0 , 0 , 0 ) 


1 u ^2 wu 3 

for any formal power series /(ci, 02 , 03 ). Substituting the dual polynomials 


ejdjaiK 


from Table |3| and the shorthands from (|21|) into (|20j) and applying c\' cq 2 C 3 3 


with i\+ 12 + i 3 = 5, we recover the Bernstein polynomials 

39 


= V WiQitf lVl + P2V2 + M 

n'*2'*3' “ 

z=i 


r n„i2„i3 
C 1 c 2 c 3 


Vi. 


Thus one sees immediately from the monomials in the dual polynomials 
which simplex splines appear in the above linear combination. For example, 
the Bernstein polynomial corresponds to the lattice vector { 11 , 12 , 13 ) = 
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(5,0,0), and 



Quasi-interpolation. For each basis B = {Si,..., S 39 } with dual points p* r , 
consider the map Q : C°( A) — > S|(A) defined by Q(f) = Yl^=i Li(f)S% and 

= E f 

k= 1 l<ri<--<r*,<5 

Note that this is an affine combination of function values of /, i.e., 



implying that Q reproduces constants. Moreover, using the Marsden identity 
it is easily checked |16:| that Q reproduces polynomials up to degree 5, i.e., 
Q(B* k ) = Bfj k , whenever i + j + k = 5. Also, Q is bounded independently 
of the geometry of A, since, using that B forms a partition of unity, 

275 

IIW)IIloo(A) < max \Li(f)\ < — ||/|| ioo(A) . 

Therefore, by a standard argument, Q is a quasi-interpolant that approxi¬ 
mates locally with order 6 smooth functions whose first six derivatives are 
in Lqo (A). Note that Q does not reproduce all splines in (A). 

Lqo stability and distance to the control points. The next theorem shows 
that each basis is stable in the L 00 norm with a condition number bounded 
independent of the geometry of A. 


+ ''' + Pi,r k ^ 


Theorem 8. Let B * = {Si, .. ., S 39 } be one of the bases B a , ■ ■ ■ ,Bf, and 
f = S T c e S| (A) with S T = (Si, ..., S 39 ) and c = (ci,..., c^) T . Then 
there is a constant > 0 independent of the geometry of A, such that 

(22) K~ 


00 — 


TA) ^ 


Proof. Let £ 1 ,..., £ 39 be the domain points of £>*. A calculation shows that 
the collocation matrix M* = { m ijffj = n with rriij = Sj(£, : ) is nonsingular, 
and its elements are rational numbers independent of the geometry of A. 

Using the Lagrange interpolant, the coefficients take the form c = Mp 1 /, 
where / = (/(£i), ■ • •, /(£ 3 g)) T - Hence, since S forms a partition of unity, 
(22) holds with Ji* = ||M“ 1 || 0 o. □ 


Note that AT is an upper bound for the condition number of the basis £?*, 
and is in fact the infinity norm condition number of the matrix M*, because 
11M* ||oc = 1- The smallest constant H-M^Hoo is obtained for B c , in which 


case 


K c = || M t 


-ii 


60866923187443943219194678615331 

836197581250152380489105335680 


72.7901. 


Hence, there is a well-conditioned Lagrange interpolant at the domain points 
of the basis B c . 

We can now bound the distance between the Bezier ordinates and the 
values of a spline at the corresponding domain points. 
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Af 

= B[0 6 0.5 1 ] 

A 4 

= A[0 5 0.5 1 ] 

A 3 := A[0 4 0.5 1 ] 

A 2 := A[0 3 0.5 1 ] 

m. 

= A[0 5 0.5 2 ] 

Bl 

= A[0 4 0.5 2 ] 

A 3 := A[0 3 0.5 2 ] 

A 2 := A[0 2 0.5 2 ] 


= A[0 4 0.5 2 l 1 ] 

Bl 

= A[0 3 0.5 2 1 1 ] 

A| := A[0 2 0.5 2 l 1 ] 
A 3 := A[0 4 0.5 2 l 2 ] 

A| := A[0 4 0.5 2 l 1 ] 

Bl 

= A[0 3 0.5 2 l 2 ] 

si 

= A[0 2 0.5 2 1 2 ] 

A 2 := A[0.5 2 l 2 ] 

Bl 

= A[0 2 0.5 2 l 3 ] 

Bl 

= A[0 4 0.5 2 1 3 ] 

A 3 := A[0.5 2 l 3 ] 

Af := A[0.5 4 1 3 ] 

Bt 

= B[0 1 0.5 2 l 4 ] 

Bl 

= A[0.5 2 l 4 ] 

A 3 := A[0.5 4 1 4 ] 


B 7 

= A[0.5 2 l 5 ] 

Bl 

= A[0.5 4 1 5 ] 



Bl 

= A[0.5 4 1 6 ] 






Table 4. Shorthands for univariate B-splines of degree d on 
the open knot multiset {0 d+1 0.5 2 l d+1 }, for d = 5,4, 3, 2. 


Corollary 1. Let h be the longest edge in A, and let f = S T c with Hessian 
matrix H and values f = (/(£i),..., /(£ 39 )) T . Then 

||c- f ||oo < 2K*h 2 max ||i?(*)|| 0 o- 
xeA 

Proof. Consider the first-order Taylor expansion of / at xq, 

f{x) - f(x 0 ) = (x- x 0 )Vf(x 0 ) + g(x). 

As the error term g G <Sf (A), it takes the form g = Y!a=\ hSi, with 

Ci - f(x o) = (£j - x 0 )Vf(x o) + bi. 

In particular for xq = $,, t . we obtain 

I ^ - /(£J| = \bf\ < AT max \g{x)\ 

*eA 

< * A* max | (a - i i ) T H{y)(x - , 

from which the theorem follows. □ 

6.2. The basis B c . While the remainder of the paper can be carried out for 
all six bases, we now restrict our discussion to the basis B c = {WiQi }?£ 1 for 
several reasons. First of all, the condition number K c is smallest. Secondly, 
this basis has the most localized support, because it contains the splines 



as opposed to splines with full support. Finally, for k = 0,1,2 and any 
direction u not parallel to the edge e of A, the number of additional splines 
Qi for which D^Qi\ e is nonzero corresponds to the dimension of the space of 
univariate splines of degree 5 — k on the knot multiset {o rf + 1_fc 0.5 2 i d + 1 ~ fc }. 
This allows for relatively pretty smoothness conditions analogous to the 
Bernstein-Bezier case. 

Derivatives on the boundary. For d = 5,4, 3,2 and i = 1,..., d + 3, let 
Bf be the univariate B-spline in Table [dj Let (ai, « 2 , a 3 ) be directional 
coordinates of a vector u with respect to the triangle [i>i, V 2 , 1 * 3 ]• Denote by 
| e the substitution of x by (1 — t)v\ + tv 2 . 

In Figure [5] and Tables !U we order the simplex splines in B c by the 
number of knots outside of e = \v\,v^. Applying ([6]), ([7]), and © we can 










1 


4Rf 

80464 

16 ot \ B \ 

2 

A- 

4 Bl 

8(a 2 6f + a \ 6 2 ) 

320402-64 + 16afR 2 

3 

A 

2 B $ 

4o2 -6 2 + 2(2ai + 02)63 

802-64 + 4 o 2(4 oi + 02)-6 2 + 2(2ai + a 2 ) 2 6f 

4 

A 

2 B \ 

2 c \ L 2 B 3 + 2(2«i + 02)64 

do^-62 + 4 o 2(2 oi + 02)63 + 2(2ai + 02 )^6 3 

5 


2 Bl 

2(ai + 2 a 2 )6 4 + 2ai 6 4 

2(oi + 202) 2 -B| + 4oi(ai + 2a2 )-B| + 4Q4-B5 

6 

A 

2 Bl 

2(ai + 202)65 + 40465 

2(oi + 2o2) 2 6 3 + 4oi (01 + 4a 2 )6f + 8af6| 

7 

A 

4 6f 

8(a 2 6| + 0:1-67) 

16a2-6| + 32aia2-6| 

8 

A 

4R| 

80:2-67 

16o 2 -6| 

9 

A 

0 

8036 4 

320403-64 

10 

A 

0 

o 3 (2R 2 4 + 6 3 4 ) 

802O364 + 03(304 + o 2 )(26 3 + -63) 

11 

A 

0 

2o3-B| 

8020363 + 203(304 + 02)63 

12 

A 

0 

4 o 3 6 4 

4a 3 (o4 + 3 a 2 )6| + 4a 3 (3a4 + o 2 )6f 

13 

A 

0 

2o 3 6 4 

203(04 + 3 a 2 )6| + 804036^ 

14 

A 

0 

o 3 (-6| + 26|) 

03(04 T 3o 2 )(6| + 26|) + 8040365 

15 

4b 

0 

8036 4 

32o2a 3 6| 

16 

A 

0 

0 

8a§6? 

17 

A 

0 

0 

2a 2 (26| + 6|) 

18 

A 

0 

0 

4a|6 3 2 

19 

A 

0 

0 

4a|6| 

20 

A 

0 

0 

2a|(6| + 26 5 3 ) 

21 

A 

0 

0 

8a 2 6| 

22 

A 

0 

0 

0 

23 

A 

0 

0 

0 

24 

A 

0 

0 

0 

25 

A 

0 

0 

0 


Table 5. Restriction of the splines Q i,... , Q 25 ; and their 
directional derivatives, to the boundary edge e = 








i 

Qi 

1 n 3 r).\ 

3-4-5 (a.\,ot.2,OL2,) Ko " l ' e 

1 

A 

32a\ B\ 

2 

* 

96a 2 a 2 Bf + 32 ofT?! 

3 

A 

80^(601 + a 2 )Bf + 4 o 2 ( 12 a! + 60102 + o!)7?| + 2(2oi + a 2 ) 3 B 2 

4 

A 

802-61 + 8 a 2 ( 3 ai + 02 )i?| + 602(201 + a 2 ) 2 B^ + 4(2oi + a 2 ) 3 B 2 

5 


4(oi + 2a 2 ) 3 B 2 + 601(01 + 2a 2 ) 2 B 2 + 80 2 (oi + 3a 2 )B 2 + 8a 3 B 2 

6 

A 

2(oi + 2a 2 ) 3 B | + 4oi (of + 60102 + 12a|)R| + 8 af (ai + 6a 2 )B 2 

7 

4b 

32a^B 2 + 96ai02-B,5 

8 

4b 

32a^B 2 

9 

* 

96a 2 a3Bf 

10 

A 

36aia 2 a3B 2 + 03(7a 2 + 5ai02 + a 2 )(2B2 + B 2 ) 

11 

A 

24o|a3 B\ + 36oi0203R| + 2os(7af + 5 qi02 + a 2 )B 2 

12 

8o 3 (o 2 

+ 50102 + 7a.2)B 2 + 803 ( 20 ^ + 7qi 02 + 2a 2 )7 ? 2 + 8a3(7o 2 + ha\a 2 + a 2 )B 2 

13 

A 

203 ( 0 ^ + 5aia2 + 7a 2 )B 2 + 360102037 ?! + 24afa37?! 

14 

A 

03(01 + 5oio 2 + 7o 2 ) (7?! + 2 7? 2 ) + 36oia2a 3 S| 

15 

A 

9602037 ?! 

16 

A 

8a|(5ai — a 2 ) B \ 

17 

A 

2402O3R1 + 2o 2 (5oi + 2o2)(27?| + 7?|) 

18 

A 

803(01 + 4 a 2 ) B 2 + 4o|(4ai + 02)7?! 

19 

A 

403(01 + 402) 7?! + 803(401 + 02)7?! 

20 

A 

2a|(5a 2 + 2ai)(7?| + 2B|) + 24 oio|r! 

21 

A 

8o|(5o 2 — oi)7?! 

22 

A 

8a|7?i 2 

23 

A 

4a|(27?2 2 + 7?3 2 ) 

24 

A 

4a|(7? 3 2 + 27? 2 ) 

25 

A 

8a|7? 2 


Table 6 . Restriction of the third-order directional deriva¬ 
tives of the splines Q i,... ,Q 25 to the boundary edge e = 
[vi,v 2 ]. 
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Figure 5 . The domain points of the bases B c and B c , or¬ 
dered by a decreasing number of knots on [111,^2], he., 1 
knot for Qi, ■ ■ ■ ,Qs, 2 knots for Q 9,..., Q 15, 3 knots for 
Q i 6 , • • •, Q 21, 4 knots for Q 22, • • •, Q25, and more than 4 knots 
for Q 26, • • • , Q 39 - 


express, for any simplex spline Q*, the restricted derivative D^Qi\ e of order 
k = 0 , 1 , 2 , 3 as a linear combination of j = 1 ,..., 8 — k. These linear 

combinations are listed in Tables [ 5 j [6] for Q 1,..., Q 25 and, by © , are zero 
for the remaining simplex splines Q 26, ■ • •, Q 39- 


Example 4 . We derive the first two entries in the third row in Table [ 5 ] 
By ©, 

area(A) 5 = 5 

area([ai,t>2,«6]) 3 3 ' 

Since u = 2 ceii>i + 20:2^4 + 2a2VQ, differentiation and knot insertion yields 




= 2ai 


+ 2a 2 


+ 2a3•0 



+ OL2 

e 

(<^D 

(2ai + a 2 ) 

+ «2 t 
e 


— 2(2«i + Oi2)B^ T 4C12-B2. 



In the next section we apply Tables [ 5 ] and [6] to derive smoothness condi¬ 
tions for splines on adjacent 12 -splits expressed in the basis B c \ see Figure [ 5 ] 
Analogously, these tables are useful for deriving smoothness conditions with 
tensor-product, Bezier, and more exotic patches. 


Smoothness conditions. Let A := [v\,V2, 113] and A := [i>i, 112, ^3] be trian¬ 
gles sharing the edge e := [vi,V2 ]■ Figure [ 5 ] shows the domain points of the 
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bases B c = {wiQifi and B c = {wiQi}i of <S 3 (A) and <Sg(A). Let 

39 39 

(23) f(v) := ^2CiWiQi(v ), v € A, f(v) := 'Y^c i w i Q i (y), v € A 

i —1 i=l 

be splines defined on these triangles. Imposing a smooth join of / and / 
along e translates into linear relations among the Bezier ordinates c,; and Cj. 


Theorem 9. Let (/5i, /? 2 , /%) &e the barycentric coordinates ofv 3 with respect 
to the triangle A. Then f and f meet with 
C° smoothness if and only if Ci = Ci, for i = 1,..., 8; 

C 1 smoothness if and only if in addition 


PlC7 + /3 2 C 8 + /?3Cl5; 

Pi Cq + P2C7 + /? 3 Ci 4 , 

/?lC5 + /3 2 (2C6 — C7) + /?3Ci3, 


2/?2/?3ClO + /?|ci6, 


C 9 = PlCl + P 2 C 2 + P 3 C 9 , C 15 

CIO = PlC2 + /32C3 + /?3ClO) C 14 

C 11 = /3i(2c 3 - C 2 ) + P 2 C 4 + focn, C 13 
C 12 = + /3 2 ^^ + P 3 C 12 ; 

C 2 smoothness if and only if in addition 
ci 6 = 0\C\ + 2 P 1 P 2 C 2 + /3f C 3 + 2/3i/5 3 c 9 
C17 = Pfc 2 + /3 2 2 C4 + /3|ci7 + 2/3 1 /3 2 3 « + 2p l p 3 ^f^ + 2(32P3 CW+2 2 11 ~ CS ’ 
CIS = /3 p C 3 +3 C4 ~ C2 + + /3|ci 8 + 2/3i/3 2 C2 ~ 2c3 + 6c4+C5 

+2/3i/ 3 3 C2 ~ 2c3+2c4 ~ 6 C5+ ^ c11+3ci2 + 2/j 2 /3 3 9ci2 ~ 2c5 ~ cl1 , 
c 19 = /3 2 ?^ + / 3 2 2C 5 + 2C6 - C7 + /3 2 ci9 + 2/3 1 /3 2 C4+6c5 ~ 2c6+C7 
+2/j 2 /3 3 C7 ~ 2c6+2c5 ~ 6 C4+3ci3+3ci2 + 2/3 i/ 3 3 9ci2 ~ 2c4 ~ C13 , 

C 20 = /?i 2 c 5 + /? 2 c 7 + /3 2 c 2 o + 2p 1 p 2 ^=^ + 2/ W 14+2 2 13 ~ C6 + 2/3 2( 5 3 ^4p^ ; 
C 21 = /3 2 c 6 + 2/?i/ 3 2 C7 + /5|c s + 2/3i/33Ci4 + 2/3 2 /53 Cis + /3|c 2 i; 

(7 3 smoothness if and only if in addition 

C 22 = / 3 ?C 1 + 3 Pl& 2 * 2 *L + 3/3i/3 2 ^=^ + p3 C4 + 3 ^ do-ca+Scn 
+3/j 2 /3g Cl(1 ~ Cl 3 6+3ci7 + /3|c 2 2 + 3/jg/3i 5ci6 3 ~ 2c9 + 

+6/3i ^3 5 ci0 ~t~ m ’ 

c 23 = /3f^a + 3 /3 2 / 3 2=2 C2±8c3±3c 4 + 2 S 2 = SLt*£4±a + 

| 3 ^ 2 ^ ciq+15ci2—C 3 — 2 c 4 —4c 5 _|_ 3 p^02 — 6 ci 2 + 2 ci 7 + 12 ci 8 —C 4 + 2 C 5 
| 3 ^ 4 ^ cio+3cn—3ci2+5ci7+3ci8+C2 — 2 C 3 +C 5 _|_ g^g2 8 ciq+3ch — 2c2 


+6/3i/3 2 /?3 


3cio+6ch+3ci2+C2 —4C3+C4—C5 


+ /^3 c 23 > 


c 24 = + 3/jf/3 2 C4+8c5 ~ C6+C7 + 3/3i/3g 3c5+8 Q 6 ~ 2c7 + 

| 3^2^ 3013+8014—207 _|_ — 3 ci 2 + 3 ci 3 +ci 4 + 3 ci<)+ 5 c 20 +c 4 —2cb+C7 

| g^^2 — 6 ei 2 + 12 ci 9 + 2 c 20 + 2 c 4 —C 5 _|_ 3 ^ 2 ^ 15ci2+ci4-4c4—2 c 5 -C6 

+6/3i/3 2 ^3 3ci2+6ci3 + 3ci V C4+C5 ~ 4c6+C7 + /3|C24, 9 

C25 = /3?c 5 + 3/3 2 /3 2 ^=^ + + /3|c 8 + 3/3 2 /3 3 ^p* 

+ 3 / 3 2 / 3 g 5 C 21 o 2C15 + 3 / 3 o/ 3 i c 14+3 C20 -c 21 + 3 ^ fl 2 3ci3+ci4-c6 

+6/3i/3 2 /3 3 5ci4 ~g 7 " C15 + /? 3 c 25 ? 

0 = (3/3g/?2 - 3/3i/3|) C2-2C3+2C4-2C5+2C6-C7 
+3/3 2 /h, C ll~3ci2+Cl3+C2— 2C3+&4 
+ 3 / 3. 2 /?3 c n~^ c i 2 + ci 3 + c 7 ~ 2 c 6+ 2c 5 
I 3 ^ 4 ^ -5cn+6ci2+ci3+9ci8-9ci9-2c2+4c3-4c 4 
| 3 0 2 p2 -5ci3+6ci2+cn+9ci9-9ci8-2c7+4c6-4c 5 

H-6/3i/3 2 /?3 ~ 2cil ~ l ~ 6ci2 ~ 2ci3 ~ c2 ~*~ 2c3 ~ 2c4 ~ 2c5 ~ l ~ 2c6 ~ c7 


Proof. By the barycentric nature of the statement, we can change coordi¬ 
nates by the linear affine map that sends ip i—>• (0,0), r? 2 1 —>• (1,0), and 
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r >3 1 —> (0,1). In these coordinates, 


«3 = /5i(0, 0) + / 3 2 ( 1,0) + /5 3 (0,1) = (( 52 , ft). 

Let u ■.= v 3 — = (/3 2 , (5%). For r = 0,1, 2,3, the splines / and / meet with 

C r smoothness along e if and only if 0) = D^f(-, 0) for k = 0,..., r. 

Substituting (23), this is equivalent to 
39 39 

^CiWiD^Qii^O) = ^CiWiDuQii-, 0), k = 0 ,...,r, 


i= 1 


Z— 1 


which, using Tables [5] and [6j reduces to a sparse system 


8 —k 


rk i B j k = °> k = 0,..., r, 

3 =1 


where is a linear combination of Cj and c* with z = 1,... ,rik+ i, where 
ni = 8 , 77,2 = 15,773 = 21, and 774 = 25. This system holds identically if 
and only if rkj = 0 for j = 1,..., 8 — k and k = 0,..., r. Let no = 0. 
For k = 0,1, 2,3, one solves for Cn k + 1 ,..., c nk+1 , each time eliminating the 
Bezier ordinates Ci that were previously obtained, resulting in the smooth¬ 
ness relations of the Theorem; see the worksheet for details 16 . □ 


As for the Bezier basis and the S-basis from |5|, each smoothness rela¬ 
tion also holds when replacing each Bezier ordinate by the corresponding 
domain point. The smoothness relations therefore also hold between the 
corresponding control points. 

The final C 3 smoothness condition only involves the Bezier ordinates in 
a single triangle and the barycentric coordinates of the opposing vertex in 
a neighboring triangle. It follows that C 3 smoothness using S 3 cannot be 
achieved on a general refined triangulation Tl 2 - 


Conversion to Hermite nodal basis. Let A be as in ([ 2 ]) with x v := Vi — v 
and y v := Vj — v for any vertex v 6 V with opposing edge [vi,Vj\ G £, 
and u e = Vk — m e for any edge e with opposing vertex v^- In addition to 
the basis B c = {wiQi} 3 ^, the spline space 5 3 (A) has the (Hermite) nodal 
basis A* = {A*} 3 ^ dual to A = {Aj} 3 £ 1 , i.e., Aj(A*) = 5ij. In this section 
we express A* in terms of B c . For details we refer to the worksheet [l6|. 

Write Qi = ajqAJ + • • • + a^gA^g for * = 1,..., 39, so that A j(Qi) = atj. 
Multiplying by the inverse of the matrix (a l j)i.j, we can express the nodal 
basis functions AJ,... A^g in terms oi Qi,..., Q 39 . 


Theorem 10. With v = v\, x = v 2 — v\, y = v 3 — v\, and u = v 3 — V 4 , 


£ 


* 

V 







+ 



+ 


1 

2 





+ 


16 


3 

+ 16 
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Note that the coefficients in these linear combinations are independent 
of the geometry of the triangle. The remaining nodal functions in A* are 
obtained by applying a symmetry in S 3 to the above equations. For instance, 



Example 5. Let v l := (cos(27ri/6), sin(27ri/6)), with i = 1, 2,..., 6 , be the 
vertices of a regular hexagon centred at the origin vq := (0, 0). Consider the 
triangulation T with triangles [t?o, v$, tq], [i>o, tq, tq],..., [^ 0 ,^ 5 , Do]- The 
nodal functions e * 0 on these triangles patch together to a spline in S^’ 3 (Tl 2 ), 
which is plotted in Figure [ 6 ] together with its control mesh. 


7. Final remarks 


Remark 11. On the Powell-Sabin 6 -split a simplex spline basis is not pos¬ 
sible, because it is not a complete graph. However, there exist C 2 quintic 
B-spline bases [1,24,25. 

Remark 12. The C 1 quadratics from 18 and the space from 


viewed as the cases n = 1,2 of a sequence of locally C 


2n-l 


14 can be 


globally C n 


spaces of degree 3n — 1. There is a natural generalization to general n of 
a set of nodal functionals that on a single triangle has size equal to the 
dimension ^n 2 + |n of this space. As remarked in 


14 


these do not form 
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Figure 6 . The nodal function e* 0 (left) and its control mesh 
(right) on the triangulation from Example [5j 

a basis for n > 2. However, it is plausible that one can instead construct 
simplex spline bases for higher degree and smoothness. 

Remark 13. A case-by-case analysis, and a computation in the worksheet, 
show that the C l quadratic simplex splines on the 12-split are 



The quadratic S-basis of A) comprises the first three types and was shown to 
have local linear independence. Taking any other combination will overload 
some of the triangles in the 12-split, and the S-basis is therefore the unique 
simplex spline basis with local linear independence. Moreover, it reduces to 
a B-spline basis on the boundary. 

For the quintic splines, a case-by-case analysis shows that each of the 
outer faces Ai, A 2 , •.., Ag will be covered by at least 9 of the 21 simplex 
splines that are nonzero on the boundary, and by at least 14 of the 18 
remaining simplex splines. Thus each outer face is covered by at least 23 > 
(5 + 1)(5 + 2)/2 simplex splines, showing that 5| admits no locally linearly 
independent S' 3 -invariant simplex spline basis that reduces to a B-spline 
basis on the boundary. 
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